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Abstract — In this paper, we propose a distributed throughput- 
optimal ad hoc wireless network scheduhng algorithm, which is 
motivated by the celebrated simplex algorithm for solving linear 
programming (LP) problems. The scheduler stores a sparse set 
of basic schedules, and chooses the max-weight basic schedule for 
transmission in each time slot. At the same time, the scheduler 
tries to update the set of basic schedules by searching for a new 
basic schedule in a throughput increasing direction. We show that 
both of the above procedures can be achieved in a distributed 
manner. Specifically, we propose an average consensus based link 
contending algorithm to implement the distributed max weight 
scheduling. Further, we show that the basic schedule update can 
be implemented using CSMA mechanisms, which is similar to 
the one proposed by Jiang et al. |1|. Compared to the optimal 
distributed scheduler in 1 1 1, where schedules change in a random 
walk fashion, our algorithm has a better delay performance by 
achieving faster schedule transitions in the steady state. The 
performance of the algorithm is finally confirmed by simulation 
results. 

I. Introduction 

Scheduling in ad hoc wireless networks is, in general, an 
NP-complete problem. In order to achieve throughput opti- 
mality, one either needs to solve a problem with exponential 
complexity in each time slot (e.g., the max-weight scheduler in 
Q), or by amortizing the complexity among an exponential 
number of time slots (e.g., the random "pick-and-compare" 
scheduler in 1,3J), so that the number of operations required in 
each time slot is constant, at the expense of large (exponential 
in the network size) queue lengths in the worst case [4 |. 

In addition, the requirement that the scheduler be imple- 
mented in a distributed manner makes the scheduling problem 
more difficult to solve. In the literature, distributed scheduling 
algorithms are often designed to be sub-optimal (e.g., the 
greedy schedulers in JS], |l6l, Q), so that minimal coordi- 
nations among links are needed during scheduling. Recently, 
Jiang et al. [1] showed a surprising result that, distributed 
throughput-optimal scheduling can be achieved by cleverly 
adjusting the CSMA contending parameters, which is based 
on the theory of Markov Chain Monte Carlo (MCMC). A 
throughput optimal scheduler is one that can achieve the 
rate stability for any arrival process whose average rate can 
be stabilized by some scheduler. This result has spurred 
interest among researchers in searching for other distributed 
throughput-optimal scheduling algorithms JU, ||9], ifTO l. How- 



ever, since these current scheduling algorithms are all based 
on the MCMC, they may suffer from large delays, even if they 
could begin in steady state (i.e., with optimal CSMA parame- 
ters chosen initially). This is due to the random walk behavior 
of the embedded time reversible Markov chain, which requires 
a long time, in general, for the schedule to change significantly. 
To illustrate this point, consider the star-shaped interference 
graph in Fig. [T] (a), where each node represents a link in 
the wireless network, and each edge represents a transmission 
conflict. Thus, the transmission schedules must be independent 
sets of the graph, e.g., {0}, and {1,2,3,4,5,6}. Due to the 
star-shaped topology, these two independent sets can also 
be viewed as the two "modes" of this network. Now, the 
distributed CSMA scheduling implies that the scheduled inde- 
pendent sets form a random walk on the family of independent 
sets, where transitions are only allowed between independent 
sets which differ by one link. For example, the schedule can 
change from {1} to {1,2}, but not from {1} to {0}, since 
that requires two simultaneous changes: shutting down link 1 
and activating link 0. Suppose the current schedule is at mode 
{1, 2, 3, 4, 5, 6}. In order to switch to the other mode {0}, all 
of the links {1, 2, 3, 4, 5, 6} need to stop transmission first, so 
that link is able to contend for transmission. Intuitively, this 
takes a long time, in particular, if the other links are unwilling 
to stop transmission, due to their own traffic demands. Thus, 
the distributed CSMA scheduling may incur a large delay by 
spending a long time in each mode before transition, since the 
random walk based design requires the scheduler to go through 
the intermediate (sub-optimal) schedules before switching to 
another mode, which happens with low probability. For more 
details about the steady state delay, please see the queueing 
simulation result in Fig. [3] which is explained in detail in 
Section |IV] 

Realizing the limitations of these distributed CSMA 
scheduling algorithms, we try to improve the delay perfor- 
mance by avoiding the random-walk based transitions between 
schedules. Instead of allocating positive probabilities to all 
the schedules, as required by MCMC, we require that the 
scheduler choose the transmitting schedules among a sparse 
set of explicitly enumerated "modes" (basic schedules), so 
that the schedule transitions are faster, since there is no need 
for the scheduler to enter the sub-optimal "transition modes". 
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Fig. 1. (a) is a star-shaped interference graph of 7 links, and (b) is a ring- 
shaped interference graph of 6 links. 

Further, since the number of basic schedules is small, optimal 
scheduling, from among these schedules, can be achieved 
efficiently in a distributed fashion. This is done by solving 
a simple max-weight independent set problem, which can be 
achieved by a certain consensus algorithm. 

The existence of such a sparse optimal solution is guar- 
anteed by the well-known Caratheodory theorem llTZI . which 
states that, at most rt + 1 schedules are needed to achieve 
the optimal scheduling in a network with n links. Further, 
such a solution can be efficiently computed by the celebrated 
simplex algorithm fTTl, which can be viewed as a steepest 
descent algorithm along the edges of the feasible region (a 
poly tope for LP). Inspired by the simplex algorithm, we 
propose a distributed algorithm to search for the optimal basic 
schedules, which can be implemented using distributed CSMA 
mechanisms similar to |[T|. Thus, by combining the consensus 
based scheduling and CSMA based basic schedule updates, we 
show that throughput optimal scheduling can be achieved in a 
distributed fashion. Finally, the delay improvement is verified 
in the simulation section. 

The organization of the rest of the paper is as follows: In 
Section we introduce the system model, and in Section |III] 
we propose and analyze the scheduling algorithm. Section |IV] 
shows the simulation results, and finally Section |V] concludes 
this paper 

II. System Model 

In this section we introduce the system model, which is 
standard in the literature. We first introduce the network model. 

A. Network Model 

We are interested in the scheduling problem at the MAC 
layer of a wireless network, where the topology of the network 
is described by an interference graph Qj — (V/, Si), where V/ 
is the set of n links, and £j is the set of pairwise interference 
constraints, i.e., link i and link j are not allowed to transmit 
together if (i, j) G £i, due to the strong interference that one 
link causes upon the other. We assume a slotted time system, 
and in each time slot t, the scheduled transmitting links crsch(i) 
must form an independent set in Qi. With an abuse of notation, 
we also denote crsch(^) as an n x 1 vector, where (Tsch,i{t) = 
1 if link i belongs to the transmitting independent set, and 
Csch.ilO = Otherwise. 



The queueing dynamics of the network is as follows: 

Q{t) = Q{Q) + A{t) - D{t) (1) 

where Q{t) is the queue length vector at time slot t, and A{t) 
and D{t) are the number of cumulative arrived and departed 
packets during the first t time slots, respectively. We assume 
that the packet arrival process is subject to the Strong Law of 
Large Numbers (SLLN), i.e., with probability 1 (w.p.l), we 
have limt_j.oo -A.{t) /t = a where a is the arrival rate vector 
Further, we also assume that the number of arrived packets in 
each time slot is uniformly bounded by a large constant. Note 
that this assumption is quite mild, since the packet arrivals are 
allowed to be correlated both across time slots, and different 
links. Thus, our model is well suited in analyzing typical 
wireless networks, and can be readily generalized to multi- 
hop scenarios. 

A basic requirement on the scheduler is the rate stability, 
i.e., for each link i, we have limj^oo = a^, w.p.l, 
so that an average throughput of a can be achieved by 
the scheduler. Thus, a throughput optimal scheduler should 
achieve rate stability for any a which can be stabilized by 
some scheduler. In the seminal paper O, Tassiulas et al. 
showed that the optimal stability region is ^ = Co(7W), 
where Co( ) denotes the convex hull, and M is the family 
of all independent sets. Further, they propose the following 
optimal max-weight scheduling algorithm: 

crsch(i) = arg max cr'^Q{t) (2) 

However, such an algorithm requires solving an NP-complete 
problem in every time slot, and is very difficult to implement 
in a distributed manner In the next subsection we introduce 
the recently developed throughput optimal distributed CSMA 
scheduling algorithm. 

B. Throughput Optimal Distributed Scheduling 

Before introducing the distributed CSMA scheduling algo- 
rithm, we need to have an optimization interpretation of the 
scheduling problem. We can formulate the scheduling problem 
as the following feasibility problem: 

SCH: minimize^; (3) 
subject to Mx = a (4) 
l^a;==l,a;^0 (5) 

where a is the arrival rate vector, M is the matrix whose 
columns are all the independent sets, and x is the scheduling 
variable, such that x^r represents the asymptotic time fraction 
that independent set cr (which is a column of the matrix M) is 
chosen by the scheduler. Thus, x naturally lives in the simplex, 
as described by (|5]l, and (|4]l is essentially the rate stability 
constraint. 

Note that SCH has many solutions, due to the key structure 
in the highly under-determined system in (|5]l that, any subset 
of an independent set is also an independent set. Thus, based 
on a solution for any strictly feasible arrival rate vector a, 
one can easily construct a set of feasible solutions. However, 



distributed implementation of any solution is a very challeng- 
ing problem. In fT\, Jiang et al. obtains a solution in the 
exponential family by transforming SCH into the following 
max -entropy problem ME: 

ME: maximizea- H(x) = — Xa-logXa- 

rreM 

subject to and (|5|) 

It can be shown that the solution x* has the following form: 

X* = exp(0*'^cr - A{e*)) (6) 

where 9* corresponds to the Lagrange multiplier associated 
with the rate constraints in dUi, and 

A{e*) = log J2 exp{e*^a) (7) 

rreM 

is the standard log-partition function. Following the literature 
of Gibbs sampling (an example of MCMC), it was shown in 
lUl that the optimal allocations x* can be implemented by 
a time-reversible Markov chain, using CSMA mechanisms. 
Further, the optimal parameters (Lagrange multipliers) 6* can 
be obtained using a stochastic gradient algorithm. In below we 
describe a discrete-time version of the distributed scheduling 
algorithm, which is also closely related to the algorithm in |8|. 



Algorithm 1 CSMA(0) 

In each time slot t, do the following: 

Randomly generate an independent set cr(t) from a certain 

distribution, independently across time slots. 

for each i e cr{t) do 

Pi = exp(6'i)/(l + exp(6'i)); 

if no neighbor of i is in <T(.sma(i ^1) then 

f,\_( 1 with probability pi 

<7csm,At) - I else 

end if 
end for 

Any other link i not in cr(i) set crc.sma,i(i) = crcsma,i{t — 1). 



In the above algorithm, carrier sensing is used in two phases: 
1) generation of the independent set cr(t) (according to certain 
protocol 1 8 1) and 2) detection of whether there is a transmitting 
neighbor of a link i G crit) at time slot t ~ 1. Further, 
note that the key part of the algorithm, phase 2, is fully 
distributed, with no explicit message exchange among links. 
The following proposition shows that the above algorithm 
achieves the distribution in (|6]l asymptotically. 

Proposition 1: crcsma(i) in CSMA(0*) form a time- 
reversible Markov chain, with the steady state distribution in 
the form of (|6]l. 

Proof: The claim is proved by checking balance equa- 
tions. The proof is essentially the same as the one in lO, 
which we omit due to space limitation. ■ 

As discussed in Section IH however, the above scheduling 
algorithm suffers from long delay, due to the random walk 



behavior associated with the time-reversible Markov chain. 
In the following section we try to solve this problem by 
introducing the simplex scheduling algorithm. 

III. Simplex Scheduling 

In this section we propose optimal distributed scheduling, 
which is based on the simplex algorithm. For a detailed 
description of the simplex algorithm for general LP problems, 
please see, for example, Since the simplex algorithm 

requires a feasible starting point, we next transform the prob- 
lem SCH into a relaxed problem REL, where a feasible initial 
vertex is easy to obtain: 

REL: minimize^. 7 (8) 
subject to Mx = (1 — 7)0 (9) 
l^x = l,a;hO,0<7<l (10) 

In above, the relaxation variable 7 can be interpreted as the 
"throughput gap", so that 7* = if and only if the arrival 
rate a can be stabilized by the scheduler. In order to illustrate 
the main ideas of the distributed scheduling algorithm, we 
first introduce a hypothesized centralized simplex algorithm 
to solve REL. 

A. Centralized Simplex Scheduling 

The centralized simplex algorithm SIM is shown in Algo- 
rithm 12] which is essentially an application of the general 
simplex algorithm to the specific scheduling problem REL. 



Algorithm 2 SIM 
Initialization: 

B = I,x*g= min( jT^, l)a,7* = max(l - jr^,0)- 
while 7* > do 

1) Simplex Search: Compute the moving direction 

Cnew = argmax l^B^^cr (11) 

rreM 

2) Scheduling: Compute the new vertex (a;^,a;*g^) and 
the throughput gap 7* by solving 

MOV: min 7 

subject to Bxb + CTnewXnew = (1 ^ T)^ 

1 X Q -\- Xj]e\v 1 

xb h 0,a;new > 0,0 < 7 < 1 

3) Update: Let be cr be a column in B such that x* = 0. 
Replace cr with cr„ew, and relabel the variables. 

end while 

return {B,Xg,j*) 



The above simplex algorithm solves the problem REL by 
moving along adjacent vertices of the feasible region in a cost 
reducing direction. In order to understand its behavior, we first 
need to identify a vertex. According to the equality constraints 
in REL, a vertex (or a basic solution in the simplex algorithm 



terminology) {xb,j) can be uniquely determined by a basis 
matrix as follows: 



B a 
1^ 



Xb 

1 



(12) 



where B is an n x n (invertible) sub-matrix of M, which 
represents n independent sets, and a;^ is an n x 1 sub-vector 
of X, which are the allocated time fractions of the independent 
sets in B. It is easy to verify that the initial vertex as specified 
in the initialization phase of SIM is feasible. 

We next show that SIM successfully moves to an adjacent 
vertex following a cost reducing direction. An edge is repre- 
sented by a new column vector (^) , where cr is an independent 
set, and moving along the edge is equivalent to increasing 
the coefficient associated with (^), with the constraint that 
the equalities in REL still hold. Thus, the changes of the 
optimization variables must stay in the null space of the matrix 
( M a> 



1^ 



i.e., we increase the coefficient of ('^) by a unit, 
the changes in the existing variables are 



B 

T 







(13) 



fAxB 
V OJ \ A7 

Using the block matrix inversion formula, we obtain the rate 
of change in the cost: 



A7 = (1 -7)(1 - l^B-^cr) 



(14) 



Thus, we need to find a proper independent set cr, so that the 
cost change A7 < 0. We next show that the cr* obtained in 
(fTTT l is a proper cost-reducing direction. 

Proposition 2: The cost change for the independent set cr* 
in (fTTT l satisfies A7 < 0, and the inequality is strict if 7 > 0. 

Proof: From (fT2l i we have xb = (1 — ^)B~^a. After 
multiplying both sides with 1^ and noting that I^xb = 1, 
we have 

(l-7)(l^B-ia) = 1 (15) 
Noting that a S Co(A^) and (fTTT l. we have 



(16) 



and therefore 



A7 = (1 -7)(1 - l^B-V*) 
< (1 -7)(1 - l^B^^a) = 



-7 



from which the claim follows. ■ 
Given the new direction specified by cr*, according to the 
standard simplex algorithm, SIM then moves along the edge 
as specified by {^-^ ), until it reaches a new vertex, where 
some coefficient of the independent set in B first becomes 
zero. Then SIM replaces that column with cr*, and relabel the 
variables if necessary. We next show that, this movement to a 
better vertex is achieved by solving the problem MOV. 

Proposition 3: For the solution {x^,x*g^,'y*) to the prob- 
lem MOV, we have x*^^ > 0, and there is one column cr in 
B such that a;* = 0. 

Proof: Due to space limitation, we only describe the 
intuition behind the proof. Suppose x*^^ — 0, then the solution 



to MOV is at the same vertex associated with the old matrix 
B, which contradicts the fact that cr^ew is a cost-reducing 
direction, according to Proposition |2] The claim that some cr 
in B has x* = follows from the fact that MOV has bounded 
optimum. ■ 

Having shown that the centralized algorithm SIM is the 
same as the general simplex algorithm for LP problems, we 
have the following conclusion: 

Theorem 1: If a G Co(A^), SIM will return a solution 
[B ,x*Q^^*) such that 7* ~ 0, and (a;Jj,7*) solves REL. 

In the next subsection we show that this algorithm can be 
efficiently implemented in a distributed manner, using CSMA 
and average consensus algorithms. 

B. Distributed Simplex Scheduling 

In this subsection we propose the distributed scheduling 
problem, and prove its throughput optimality, by showing that 
it is a distributed implementation of the centralized simplex 
scheduling algorithm in the last subsection. We first illustrate 
the full scheduling algorithm DIS-SCH in Algorithmic] 

Algorithm 3 DIS-SCH 

Initialization: Set the initial parameters as -6(0) = /, 

0(0) = 0, 7(0) = 1, (Tesma(O) = 0, and <T„ew(0) - 0. 

In each time slot t, do the following: 

1) CSMA: update (T^^^.^{t) by running CSMA(a0(t)), 
where a is a sufficiently large constant. 

2) Scheduling: for each link i, compute 

^.W(i) 



0'sch,i(i) 



arg max 
<Te_Bu{<T„, 



Wo 



.} 



(17) 



Link i transmits if CTsch,j:(^) = 1- In above, Wa-{t) — 
6{t)'^(7 is the weight of independent set cr, and 
w^' (t) is link i's local copy, which is updated by the 
consensus algorithm. 
3) Consensus: update the parameters 

eit) = 0(t)+e((l-7(i))a(<)-rT,eh(t))(18) 



7(i) 
a{t) 



[7(t) + e(0(t)^a(t)-l)][o,i] (19) 



A{t)/t 



(20) 



where [■]io,i] means projecting onto the interval [0, 1], 
and e is a constant step size. Run an average consensus 
algorithm over the quantities {wa-{t)}rT^Bu{rT } and 
lit). 

4) Update: If certain convergence conditions are satis- 
fied, replace the column in B with the smallest weight 
by Cnew(i), and load the CSMA independent set by 
letting cr„e„(t) = Crcsma(i)- 

Overall DIS-SCH is very similar to SIM, with the following 
two important changes: 1) the simplex search phase in SIM 
is replaced with the CSMA phase, and that 2) the centralized 
moving procedure MOV is replaced with a simple max-weight 
scheduling in ( fTTT ). which is then followed by parameter 
updates and average consensus algorithm. It may appear that 



([TT) is similar to the max-weight scheduler in (|2), with 
the difference that the queue lengths Q{t) are replaced by 
the "virtual queue lengths" 6{t). However, the max-weight 
independent set in ( [TtT i is chosen from n + 1 independent 
sets, instead of the family of all independent sets, which 
has exponential size. Thus, ( fTTl ) can be solved efficiently 
(in linear time), whereas the max-weight scheduler in (|2) is 
NP-complete, in general. In the following we elaborate on 
the change in Step 2, by showing that the scheduling and 
consensus phases in DIS-SCH is equivalent to solving the 
problem MOV. 

Proposition 4: If there is no change in the matrix B, the 
parameters d updated in ( fTSl l will converge to the optimal 
dual variable for MOV, 7 will converge to the optimal cost 
7*, and the average time fractions Xnew) will converge to 
the optimal primal variables (a?^,a;*g^) for MOV. 

Proof: Due to space limits, we only describe the intuition 
behind the proof. We first form the Lagrangian of MOV 

minimize3;3,a;„e„,7 7 + ^^((1 - 7)a - Bxb - cTnewa^new) 
subject to l^ccs + x^ew = 1, ^0, a;new > 

Note that this is a LP over a simplex, and therefore, the 
solution is obtained at a vertex. Specifically, we can choose 
the vertex Xo-,^,, — 1 (which is equivalent to scheduling 
independent set cTsch), where satisfies 

Csch G arg max a (21) 

Further, under mild assumptions about the consensus algo- 
rithm, it can be shown that the two schedules in ( [TtT i and (ISTT i 
are the same with high probability. Finally, it can be shown that 
the updates in (fTsT l and (fT9] l correspond to the standard primal- 
dual algorithms for solving convex optimization problems. 
Thus, the variables will converge to the optimal. ■ 
We next elaborate on the change in Step 1, by showing 
the equivalence between the CSMA phase in DIS-SCH and 
the simplex search phase in SIM. We have the following 
proposition. 

Proposition 5: With sufficiently large a and sufficiently 
long time, <J^^m3.(f) is the max-weight independent set in ( fTTT l. 
with high probability. 

Proof: Due to space limits, we only describe the intuition 
behind the proof. From the complementary slackness property, 
from the Lagrangian of MOV we have 

e'-^B = (1 - 7*)(0*a)l^ (22) 

Note that for notation simplicity, B has already been up- 
dated by replacing one sub-optimal column with cTnew, i e., 
B = i?new Thus, the optimal dual variables 0* satisfies 
Q* oc B^^ , and therefore, if a^^^JyC) satisfies 

cTcsm^it) e argmax0*^cr (23) 

then it is essentially the same max-weight solution of (fTTT l. 
Further, note that the CSMA based sampling phase converges 
slowly (exponential time, in general), whereas solving MOV 
are relatively fast. Thus, we can assume that the parameters 



6{t) have converged to the optimal, and finally, the claim 
follows from the steady state distribution in and the well 
known approximation that 

exp(a6>*^cr - A{a9*)) « l{^gargmax„eM (24) 

for large enough a, where is the indicator function, i.e., 

l{true} = 1 and l{ake} =0. ■ 

Having shown that the distributed algorithm DIS-SCH is 
essentially the same as the centralized simplex algorithm SIM, 
we have the following conclusion: 

Theorem 2: If a G ^ DIS-SCH will return a solution 
(i3,cc^,7*) such that 7* = 0, and (a;Jj,7*) solves REL. 

In the next section we will demonstrate the performance of 
DIS-SCH by simulation results. 

IV. Simulation Results 

We next compare the performance of the distributed simplex 
scheduling algorithm DIS-SCH with the distributed CSMA 
scheduling algorithm in Algorithm [T] using MATLAB simula- 
tion. During the simulation, we assume that the packet arrivals 
are i.i.d, with 95% of the maximum uniform arrival rate. 
The average consensus algorithm for the distributed simplex 
scheduling is implemented as follows: We assume that the 
communication graph for the average consensus is the same 
as the interference graph. In each time slot, a random maximal 
matching is first formed, and each matched pair update their 
variables by taking the average of their local copies. We also 
allow a certain time period for the average consensus algorithm 
to converge, before it is used for distributed scheduling. 

A. Star Network 

We fist consider the 7-link star-shaped interference graph 
in Fig [U (a), with the simulation result shown in Fig. |2] 
In the figure, the sample queue lengths are plotted over a 
simulation period of 2 x 10^ time slots. From the figure, one 
can observe that the network is rate stable in both cases, but 
distributed CSMA scheduling has much larger queue lengths 
(around lO"*) than simplex scheduling (several hundreds) in 
the steady state. Further, one can observe that link is the 
bottle neck Unk for CSMA scheduling, since its queue length 
is most often the largest. This is because, as discussed in 
Section H] in the steady state, distributed CSMA scheduling 
spends a considerable amount of time around each mode 
before transiting to the intermediate (suboptimal) independent 
sets. Thus, the transitions of CSMA scheduling is very slow, 
and the queue lengths are consequently large. On the other 
hand, simplex scheduling can quickly switch between the two 
optimal modes, and therefore, have smaller queue lengths in 
the steady state. 

B. Ring Network 

We next consider the 6-link ring-shaped interference graph 
in Fig [U (b). The simulation result is shown in Fig. |3] 
Similar to the star network, one can observe that both al- 
gorithm result in rate stability, but the simplex scheduling 
achieves much smaller queue lengths than the distributed 
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Fig. 2. The simulation result of a 7-star network with distributed CSMA 
scheduling and simplex scheduling. 
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Fig. 3. The simulation result of a 6-ring network with different distiibuted 
CSMA scheduling and simplex scheduling. 



CSMA scheduling. In fact, for the ring shaped network, it is 
easy to see that the optimal modes are cri = (1, 0, 1, 0, 1, 0)"^ 
and (72 = (0,1,0,1,0,1)^. In the steady state, the simplex 
scheduler can achieve low delay by quickly switching between 
these two modes, using the simple max-weight scheduler 
implemented by average consensus. On the other hand, the 
switching is much more time-consuming for the distributed 
CSMA scheduling, due to the random walk based design. 
Finally, note that the transition time and queue lengths for 
distributed CSMA scheduling in the 6-ring network is smaller 
than that for the 7-star network (both have the same arrival 
rates), due to the fact that, it is easier to switch between the 
modes in a ring-shaped topology. 



V. Conclusion 

In this paper, we proposed a distributed throughput-optimal 
scheduling algorithm for ad hoc wireless networks, which is 
motivated by the simplex algorithm for solving LP problems. 
The scheduler maintains a sparse set of basic schedules, and 
during scheduling, the basic schedule with the maximum 
weight is selected for transmission in each time slot. The set 
of basic schedules are updated according to the simplex algo- 
rithm, which can be implemented using CSMA mechanisms 
in a distributed fashion. Compared to the distributed CSMA 
based scheduling in ([T], our algorithm achieves better delay 
performance in the steady state by allowing faster transitions 
among the optimal independent sets. 
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